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Abstract 

^ \ Superconductivity (SC) may microscopically coexist with density wave (DW) when the nest- 
r^" 1 ' ing of the Fermi surface (FS) is not perfect. There are, at least, two possible microscopic 
q ■ structures of a DW state with quasi-particle states remaining on the Fermi level and leading 
^ . to the Cooper instability: (i) the soliton-wall phase and (ii) the small ungapped Fermi-surface 
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pockets. The dispersion of such quasi-particle states strongly differs from that without DW, 
and so do the properties of SC on the DW background. The upper critical field H c2 in 
such a SC state strongly increases as the system approaches the critical pressure, where 
superconductivity first appears. H c2 m ay considerably exceed its typical value without DW 
and has unusual upward curvature as function of temperature. The results obtained ex- 
plain the experimental observations in layered organic superconductors (TMTSF) 2 PF 6 and 
a-(BEDT-TTF) 2 KHg(SCN) 4 . 
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^ 1. Introduction 

The interplay of superconductivity (SC) and density-wave (DW) states appears in many 
00 ! compounds and is extensively investigated for several decades (see, e.g., reviews in Refs. 

[H,Ei]). The DW removes electrons from the Fermi level, and, usually, precludes superconductivity^ 
or strongly reduces the SC transition temperature [4]. However, in several compounds 
^ . in layered organic superconductors (TMTSF) 2 PF 6 and a-(BEDT-TTF) 2 KHg(SCN) 4 , 
^ \ the SC transition temperature T^ c on the DW background is very close to (or even exceeds) 
T<f without DW. In both these compounds SC coexists with DW in some pressure interval 
P c i < P < P c : superconductivity first appears at P = P c x, and at P, the DW phase un- 
dergoes a phase transition into the metallic state (see Fig. 7 in Ref. [5( and the schematic 
phase diagram in Fig. [I]). In this region SC state has many unusual properties, such as 
the divergence of the upper critical field H c2 at P — > P c i-@, 0] To explain this property, the 
scenario of macroscopic spatial separation of SC and DW states was proposed, where the 
size d s of SC domains depends on magnetic field and in strong magnetic field becomes much 
smaller than the SC coherence length When the width d s of a type II superconducting 

slab becomes smaller than £sc, the upper critical field H c2 is increased compared to H® 2 in 
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Figure 1: The schematic picture of the phase diagram in (TMTSF^PFg, where superconductivity coexists 
with DW in some pressure interval above P c i but below P c . DWo stands for the uniform fully gapped DW. 
DWi denotes the DW state when the imperfect nesting term 2t' b > A , so that the ungapped FS pockets or 
nonuniform DW structure appear. 



a bulk superconductor by a factor (see Eq. 12.4 of Ref. |8|) 

H c2 /H° c2 « VU£sc/d a . (1) 

For (TMTSF) 2 PF 6 and for many other DW superconductors this scenario implies the size of 
the SC domains to be of the order of the DW coherence length t^owi which means a strong 
influence of the DW order parameter on the SC properties and a microscopic coexistence of 
these two states. Such a strong spatial modulation of the DW order parameter costs energy 
W of the order of the DW energy gap A , which is much larger than the energy gain due 
to the surviving of SC state. In this case, the soliton wall scenario^, QiJ 11, 12, 13, 14] is 
more favorable, where the energy loss W due to nonuniform DW order parameter is compen- 
sated by the large kinetic energy of soliton-band quasiparticles. [HI, [l3| SC appearing in the 
soliton wall phase of spin-density wave (SDW) has the triplet SC pairing|l4j] in agreement 
with experiments in (TMTSF^PFg. 15] Another microscopic structure was proposed 14] and 
investigated 16|] recently, where the small ungapped Fermi-surface pockets in the DW state 
appear due to the imperfect nesting of Fermi surface (FS) and lead to the SC instability. 
In that scenario the divergence of upper critical field H c2 at P — ► P cl is due to the de- 
crease of the mean square velocity on the Fermi surface when the size of the new FS pockets 
decreases. [l6| In the present paper, we continue to study the SC properties on the DW back- 
ground. We calculate in the soliton-wall scenario and show how the divergence of 
at P — > P c i appears on the soliton-phase background. Within this scenario we also explain 
the upward curvature of the temperature dependence of H c2z along the z-axis perpendicular 
to the conducting layers. The possibility to explain the same feature within the open-pocket 
scenario and the limitations of the weak-coupling theory to describe this effect are briefly dis- 
cussed. We also explain the observed hysteresis in the pressure dependence of SC transition 
temperature T c (P) in both above scenarios. 
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2. The model and the DW state without superconductivity 



Below we consider the same model, as in Ref. 14j. The quasi-lD free electron dispersion 
without magnetic field has the form 

e(k) = ±v F (k x Tk F )+t ± (k ± ), (2) 

where the interchain dispersion t±(k±) is much weaker than the in-plane Fermi energy v F k F 
and given by the tight-binding model with few leading terms: 

t±(k±) = 2t b cos(k y b) + 2t' b cos(2k y b). (3) 

Here b is the lattice constants in the ^/-direction. The dispersion along the z-axis is con- 
siderably weaker than along the y-direction and is omitted. The FS consists of two warped 
sheets and possesses an approximate nesting property, e(k) ~ —e(k—Q N ), with Q N being 
the nesting vector. The nesting property leads to the formation of DW at low tempera- 
ture and is only spoiled by the second term t' b (k y ) in Eq. (j3J), which, therefore, is called 
the " antinesting" term. Increase of the latter with applied pressure leads to the transi- 
tion in the DWi state at P > P c \, where the quasi-particle states on the Fermi level first 
appear and lead to the SC instability. In the pressure interval P c \ < P < P c the new 
state develops, where the DW coexists with superconductivity at rather low temperature 
T < T C 5C , while at higher temperature Tsc < T < T^w the DW state coexists with the 



metallic phase. This coexistence takes place via the formation of small ungapped pockets 16 



or via the soliton phase |12l. Il3l . Il4j . We take the DW transition temperature to be much 
greater than the SC transition temperature, T® w 3> Tj? , which corresponds to most DW 
superconductors. For example, in (TMTSF) 2 PF 6 T^ DW « 8.5K > T$ c « I.IK, and in 
a-(BEDT-TTF) 2 KHg(SCN) 4 , T? DW « 8K > T c sc « O.IK. Therefore, we first study the 
structure of the DW state in the pressure interval P c \ < P < P c , and then consider the 
superconductivity on this background. 

In the soliton phase, which presumably establishes in the pressure interval P c \ < P < P c , 
the DW order parameter depends on the coordinate along the conducting chains: A (x) ~ 
A sn (x/£dw)i where sn(y) is the elliptic sinus function. As a result an array of soliton walls 
of width £dw get formed, where the DW order parameter changes sign. Each soliton wall 
contributes one electron-like quasiparticle per conducting chain on the Fermi level, which 
leads to the formation of the new conducting soliton band. This soliton band appears in the 
middle of the DW energy gap and has the dispersion 

E(k) = E(k x )+e + (k y ), (4) 

where the interchain part of the dispersion is given by the antinesting term in the dispersion 
(J3J): £ + (k^) = \t i (k i ) + tj_(k_|_ — QjJ] /2. The dispersion E {k x ) along the conducting chains 
was found in Ref. [l9|] (see Fig. 1 in Ref. [19(), and for qualitative analysis it can be 
approximated by 

E(k x ) fxE-an[n(\k x \ - k F ) /2k ] . (5) 
The soliton band width E_ and boundary kq in the momentum space are determined by the 



linear concentration n Q of the soliton walls 19]: kq = 7rn s /2, hv F n s = E + /K ( yj 1 — E 2 _ / E\ J , 
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where K(r) is the complete elliptic integral of the 1st kind. At n s — > E + k, Ao, and E- ~ 
4A exp ( — A / hvFn s ). Note, that the soliton phase can be energetically most favorable DW 
state at imperfect nesting. 13] Compared to the uniformly gapped DW state, the soliton 
phase gains the kinetic energy of quasiparticles in the half-filled soliton phase due to the 
term s + (k y ) in which compensates the energy loss of the non-uniform DW structure. 

The soliton wall concentration n s depends on external pressure P via the dependence on 
electron dispersion. n s (P) can be obtained by minimization of the total energy of the soliton 



phase (see Eqs. (33)-(35) of Ref . [ll| or Eqs. (9)-(12) of Ref . [l3|p. The result substantially 
depends on the harmonic content of electron dispersion t± (ki )[13j]. For typical dispersion, 
at P — > P c i the soliton band width E^/Aq — > 0, which corresponds to a phase transition 
from the uniform DW to the soliton phase. In the vicinity of P = P c i, rough estimates 
give 
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A <5 = a/Aq (2t' b - A c 



p-p 



cl- 



(6) 



3. SC properties on the soliton-phase DW background 



In the Ginzburg-Landau functional for SC state the coefficient tensor l/my before the 
gradient term can be expressed via the product (v iVj) of electron velocities Vi averaged over 
the Fermi surface [2oj] [see Eqs. (48)-(50) of Ref. jig]. Using this result and the simplified 
dispersion (jljhflSJ) in the soliton band one easily obtains 
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along the z-axis is 
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where ^ DW = 
44.2. Eq. (JE! 
size 8 = 2t' b 



hv F /7rA , n s £ DW w l/Vm (4x/2A /E_), and the constant C ls = 12vr 3 /7C (3) 



is similar to Eq. (55) of Ref. 
— Aq is replaced by the soliton band width E 



161 ] for the second scenario, where the open pocket 
At P -> P c i both 5 -> and 



EL — > 0, and in both scenarios the functions 5 (P) and E^ (P) depend on the electron dis- 
persion t± (kj_). Substituting (jSj) into dHJ, one obtains that in the soliton- wall scenario the 
slope dH z c2 jdT oc 1/y/P — P c \ and the upper critical field H c2 diverge at P — > P c i similar to 
the scenario of the ungapped FS pockets [161] . The physical reason for this divergence is the 
strong change of quasi-particle dispersion on the Fermi level when P — > P cl and the soliton 
band shrinks. This shrinking leads to the decrease in the mean square electron velocity on 
the Fermi level (or diffusion coefficient), which increases the upper critical field H C 2. 
For (TMTSF) 2 PF 6 the substitution of £L ~ ^/A^5 to (0) gives the slope 



dH*2 



7.1 



T, 



sc 



dT In (4 V / 2A^75) v 7 ^ 



Tesla 
°K 



(9) 
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and the maximum slope dH z 2 /dT pa 0.6 [Tesla/°K] in a reasonable agreement with the 
experiment (see Fig. 2 in Ref. Q). Far from P = P c i, £L ~ Ao, and Eq. (|Sj) gives 
dH z 2 /dT pa 0.25 jTesla/°K], which is again in a good agreement with experimental data in 
Fig. 2 of Ref. Q. For a-(BEDT-TTF) 2 KHg(SCN) 4 , the substitution of P_ « A to (JED, 
using BCS relation A = 1.76T° DW , gives H z c2 pa 7.7mTesla- (l - T/T C 5C ), or dH z 2 /dT pa 
77 \mTesla/° K] in a qualitative agreement with experiment 0] . 

The upward curvature of H Z 2 (T) in the scenario of soliton walls is similar to H c2 (T) 
in layered superconductors when magnetic field is parallel to the conducting layers. The 



latter case was considered in a number of theoretical papers (see, e.g., [21[ , [22] , [23J , or §16 
of 0). The ID network of soliton walls can be considered as a ID Josephson lattice, where 
the conducting layers of thickness £dw are separated by the insulating layers by thickness 
s = l/n s . At low field this Josephson lattice behaves like a 3D superconductor with critical 
field H l c2 = ejjjfc$ /27r£j (T) £ fc (T) , where $ is magnetic flux quantum and £j (T) is the SC 
correlation length in i-direction. In the vicinity of critical temperature & oc (T c sc — T) -1 / 2 . 
At high field H c2 tends to its value ([TJ in the 2D SC slab of thickness d s ~ ^dw Csc-I^l 
For (TMTSF) 2 PF 6 , £ DW /£ SC ~ 1/10, and according to Eq. CQ), the upper critical field E z c2 
on the SDW background can be enhanced by a factor H c2 /H® 2 < 30 as compared to H* 2 in 
the bulk superconductor, i.e. superconductor without the soliton structure. 

The crossover from the 3D to 2D behavior of H z 2 occurs when the soliton bandwidth be- 
comes small: E- < h 2 /m y ^g C , where £sc ~ hv ymax /7rAsc (T) is the temperature dependent 
correlation length within the conducting layer, and m y pa t' b /v yraSuX . This gives the crossover 
value E_ ~ [vrA^c* {T)] 2 /t' b . Therefore, the upper critical field increases and behaves as in 
the isolated SC slab only at P — > P c \ and only at low temperature, which means the unusual 
upward curvature of H z 2 (T) at P — ► P c \. Note, that the soliton structure is important for the 
upward curvature of only z-component H z 2 (T) of magnetic field in layered organic metals as 
(TMTSF) 2 PF 6 or a-(BEDT-TTF) 2 KHg(SCN) 4 , because it creates the layered soliton wall 
structure parallel to the magnetic field. If magnetic field lies in the x-y plane, it is already 
parallel to the conducting molecular layers, and H Z 2 {T) may have the upward curvature 
according to the Lawrence-Doniach model even without soliton walls. 

4. Discussion and summary 

The above study gives a rough estimates of the upper critical field H c2 (P, T) in the 
soliton phase. The divergence and the unusual upward curvature of H z 2 (T) at P — > P c \, 
observed in (TMTSF) 2 PF6 and a-(BEDT-TTF) 2 KHg(SCN)4, are explained. A quantitative 
study of H c2 (P, T) requires a more accurate model for the coexisting SC and DW states, 
which includes the influence of SC on the DW state. The critical fluctuations near the DW- 
metal phase transition also strongly influence the SC state via the renormalization of the e-e 
interaction, [l], 0, 25] However, the proposed rough model gives a reasonable quantitative 



agreement with experiment on the slope of H z c2 (T) at T — > T C 5C . 

One more unusual feature observed in these compounds is the hysteresis in the pressure 
dependence of the SC transition temperature (P). In the soliton scenario this hysteresis 
may come from the motion of soliton walls in order to achieve the optimal soliton-wall 
concentration n s (P) as pressure P changes. The hysteresis in T^ c (P) can also be explained 
in the open-pocket scenario if the DW wave vector shifts at P > P& with change in pressure 
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to find the optimal DW wave vector at imperfect nesting. Since the DW wave vector is 
pinned by crystal imperfections, this shift of the DW wave vector is hysteretic, which leads 
to the hysteresis in S and in the SC transition temperature T$c- 

An accurate description of the SC properties in both scenarios must go beyond the 
weak-coupling theory, because the new Fermi energy (i.e., 5 in the new FS pockets or the 
bandwidth E_ in the soliton phase) is comparable to the SC energy gap as P — > P c \ and 
is much smaller than the Debye energy. This fact, for example, may be responsible for the 
upward curvature of H* 2 (T) in the scenario of open pockets, where 

H z c2 oc l/yj{v*) {vl) oc 1/5. (10) 

When 5 < T C 5C , the small value of 5 is smeared out by high temperature. Qualitatively, 
this smearing may be included via the replacement 5 — > a/<5q + aT 2 , a ~ 1. At P — > P& 
substituting this and So oc P — P c \ to ( fTOl) . one obtains 

H z c2 oc l/^(3{P-P cl f + aT\ (11) 

which gives the upward curvature of H* 2 (T) and describes well the experimental data in 
Ref. |7| with two fitting parameters a and (3. 

The proposed approximate model of the SC state on the soliton-wall background of 
the DW at P > P c \ qualitatively explains all unusual SC properties on the DW background 
observed in the organic metals (TMTSF^PFg and a-(BEDT-TTF)2KHg(SCN)4, which gives 
a strong argument in favour of the microscopic coexistence of superconductivity and density- 
wave state in these compounds. 

The work was supported by RFBR No 06-02-16551, No 06-02-16223, by MK-4105.2007.2 
and the research program of JSC "Advanced Energetic Technologies". 
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